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ABSTRACT 


We clarify a relationship of the dynamics of a solar flare and a growing Coro¬ 
nal Mass Ejection (CME) by investigating the dynamics of magnetic helds during 
the X2.2-class flare taking place in the solar active region 11158 on 2011 February 
15, based on simulation results obtained from Inoue et ah (2014a). We found 
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that the strongly twisted lines formed throngh the tether-cntting reconnection 
in the twisted lines of a nonlinear force-free held (NLFFF) can break the force 
balance within the magnetic held, resnlting in their lannch from the solar snrface. 
We farther discover that a large-scale hux tnbe is formed dnring the ernption as 
a resnlt of the tether-cntting reconnection between the eruptive strongly twisted 
lines and these ambient weakly twisted lines. Then the newly formed large hux 
tube exceeds the critical height of the torus instability. The tether-cntting re¬ 
connection thus plays an important role in the triggering a CME. Furthermore, 
we found that the tangential helds at the solar surface illustrate diherent phases 
in the formation of the hux tube and its ascending phase over the threshold of 
the torus instability. We will discuss about these dynamics in detail. 


1. Introduction 


Solar hares and coronal mass ejections (CMEs) are powerful explosions of the solar 
coronal plasma that are widely considered to be sudden liberation of the free magnetic 


energy in the solar atmosphere (Forbes 2000 Priest & Forbes 2002). Our understanding of 


them is important so that we can not only clarify the nonlinear dynamics of the solar plasma, 
he., the store-and-release processes for magnetic energy and heilicty but also establish the 


space weather forecast. Since |Carrington 1859 discovered the solar hare, many observational, 
theoretical, and numerical studies have been conducted in order to clarify the properties of 


solar hares and their relationship with CMEs (Benz 2008; Shibata & Magara 2011). 


Carmichael 

(1964 

), 

Sturrock 

(1966 

), 

Hirayama 

(1974 

) and 

Kopp & Pneuman 

(1976) 


constructed a hare model based on magnetic reconnection, which well explained observations 
of solar hares made at multiple wavelengths. This model is named the CSHKP model and 
is today considered a standard hare model. Since this model was formed, solar physics 
satellites have provided countless data along with images of the solar hares and CMEs, 
which are further clarifying their behaviors. In particular, the Yohkoh satellite convinced us 


of a hare scenario based on the magnetic reconnection (Tsuneta et ah 1992 Masuda et ah 


1994). In addition, Yokoyama et ah (2001) discovered clear evidence of reconnection inhows 


associated with the solar hare from the observations of Yo/ifcoh/SXT and SOHO/YYI. All 
of these observations built up a foundation of the reconnection model for the solar hare 
as well as support the classical CSHKP hare model. These observations and theoretical 
models clarihed the properties of solar hares and helped us understand hare physics in a two 
dimensional regime. On the other hand, our understanding of the dynamics of solar hares is 
not fully developed in a three dimensional (3D) regime. 
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Recently, an innovative development of a computer allows us to begin closing this de¬ 
ficiency in our knowledge by performing a 3D magnetohydrodynamic (MHD) simulation to 


explore solar flare dynamics in 3D space. Antiochos et ah (1999) performed 3D simulation 


and propose a magnetic breakout model showing an onset of the CME by removing a overly 
field lines above the core field, through a reconnection. The model has been extended and 
high resolved simulations are recently performed {e.g., Lynch et al.|2008 Karpen et al.|2012). 


Amari et ah (2000) and Amari et ah (2003) successfully showed the formation of a flux tube 


and its eruption via a twist and converging motion on the photosphere from an initial state 


which starts from a potential held. Torok & Kliem (2005) conhrmed that the helical kink 


instability well explains the initial processes of a CME started with anchoring a hux tube in 
the photosphere and embedding it in the solar corona. They further pointed out that the 
decrease in the overlying held is a key factor determining whether the resulting eruption fails 


or is ejective. Forbes (1990) constructed the hux tube model in 2D space and pointed out 


the eruption via the loss of equilibrium. Following 

Forbes 

(1990 

), 

Inoue & Kusano 

(2006) 

applied a simple straight hux tube in the 3D solar corona; however. 

Inoue & Kusano 

(2006) 


indicated that the hux tube exists in an unstable state against the kink instability before 
it loses its equilibrium. They further suggested an existence of a critical height deciding 
whether the eruption is full or failed. Aulanier et al. 2010 and Fan (2010| well-reproduced 


the initiation of a CME driven by the torus instability (Kliem & Torok 


Aulanier et 


al. (2012) brought an interpretation of the CSHKP model into 3D space, explaining the 3D 


ehect, e.g., the strong-to-weak shear transition in the post-hare loop which is not seen in the 


classical CSHKP model. Kusano et al. 2012 explained that two diherent types of emerging 
small huxes play a role in the triggering mechanisms for the solar hare and CMEs. Recently, 
a more realistic situation is constructed; where the initial condition mimics a real situation 
or is applied to a potential held extrapolated from the vector held, enabling us to compare 
3D simulations with real observations {e.g., Zuccarello et al. 2012; van Driel-Gesztelyi et al. 


2014). More recently. An & Magara (2013), Magara (2013) and Leake et al. (2014) are try¬ 


ing to understand the onset mechanism from their results on the hux emergence simulation 
without constructing a model. 

These 3D simulations show complicated and specihc nonlinear dynamics only observable 
in 3D space. In addition, they interpolate well between 3D dynamics and the observations or 
models constructed in 2D space. Recently, a solar optical telescope ( Tsuneta et al.||2008 ) and 
the Helioseismic and Magnetic Imager ( Scherrer et al.| 2012 Hoeksema et al. [2014 ) on board 
the latest solar physics satellite Hinode ( Kosugi et aL||2007[ ) and Solar Dynamics observatory 
(SDO: Pesnell et aL||2012 ) have provided vector magnetic helds in high temporal and spacial 
resolutions. These vector helds allow us to reconstruct the 3D coronal magnetic held under 


a nonlinear force-free held (NLFFF) approximation (Wiegelmann & Sakurai 2012 Regnier 
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2013) which takes into account even the tangential components of magnetic helds of the 


vector magnetic helds and is much different from the potential held extrapolated only from 
the normal component of the vector magnetic helds. The NLFFF exhibits the strong sheared 
magnetic helds appearing above the polarity inversion line (PIT) prior to the solar hares and 
CMEs {e.g. Regnier & Amari 2004 or Canon & Amari 2010), as observed by ground and 


space observations. These are not shown in the potential held. On the other hand, it shows 


the relaxed magnetic held lines after the hare (Schrijver et ah 2008 


2008; 

Inoue et ah 

2012; 

He 


et ah 2014). Some papers further clear the temporal evolution of the magnetic held in the 


active region throughout the hare event ( 

Thalmann & Wiegelmann 

2008; 

Inoue et ah 

2011 


Sun et ah 2012 Inoue et ah 2013 Jiang et al. 2014 Cheng et al. 2014). However, because 


these NLFFFs are constructed in the steady sate, they never reveal the dynamics of the 
solar hares. Therefore, the MHD simulation which takes into account the NLFFF, i.e., data- 


constrained simulation, are taking a new approach (Jiang et al. 2013 Kliem et al. 2013 


Pagano et al. 2013). 


More recently, Inoue et al. (2014a) performed the MHD simulation using the NLFFF 


as an initial condition in order to elucidate the dynamics of the solar hare. This NLFFF is 
reconstructed 2 hours approximately before the X2.2-class hare in AR11158. The AR11158, 
which consists of a complicated quadrupole structure, produced several M- and X-class hares 


in 2011 February (Schrijver et al. 2011 Sun et al. 2012; Inoue et al. 2013) and Hinode and 


SDO satellites provided rich data observed in multiple wavelengths over the course of the hare 
events. Inoue et al. (2014a) showed that the dynamics obtained from our MHD simulation 


successfully reproduced observed phenomena such as the distribution of the two-ribbon hare 
and the held lines structure as seen in the extreme ultraviolet (EUV) post-hare image taken 
by Atmospheric Imaging Assembly (AIA; Lemen et al. 2012| ) which was on board SDO. 
However, the detailed dynamics of this event are not yet clear. 

In this study, we explore the detailed dynamics of the magnetic held during the X2.2 


solar hare taking place in AR 11158, based on the simulation results obtained by Inoue et al. 


(2014a). In particular, we clarify the dynamics connecting the solar hare and the CME, i.e. 
the geometry and evolution encompassing the region of the hare and CME initiation. This 
paper is constructed as follows. Our observational data and numerical method are described 
in Section 2. Our results and discussions are presented in Sections 3 and 4. Our conclusion 
is summarized in Section 5. 
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2. NLFFF Extrapolation and MHD Simulation Methods 
2.1. Numerical Method and Observations 


All of the methods for the NLFFF and MHD simulation following Inoue et ah (2014a). 
Therefore, here we briefly provide an overview of the numerical methods and observations 
presented in that work to aid in readers understanding. The NLFFF and MHD simulations 
are based upon the following equations; 


P = \B\, 


di) 1 

— = —iv ■ V)v H —J X B + u'V'^v, 
ot p 


OB 

— = Vx(nxS - rjiJ) - V0, 


( 1 ) 

( 2 ) 

( 3 ) 


J = V xB, 


( 4 ) 






( 5 ) 


where the subscript i of rj corresponds to ’NLFFF’ or ’MHD’. The length, magnetic held, 
density, velocity, time, and electric current density are normalized by L* = 216 Mm, B* = 
2500 G, p* = \B*\, VX = B* /{pop*y^^, where po is the magnetic permeability, = L*/VX, 
and J* = B*/pqL*, respectively. The non-dimensional viscosity u is set as a constant 
(1.0 X 10“^), and the coefficients c^, ci in equation (O also hx the constant values, 0.04 


and 0.1, respectively. Although the density is given as a model like Amari et ah ( |1996 ), we 
previously discussed its validity in [Inoue et ah (2014a). 

The vector magnetic held that we used as the boundary condition is observed at 
00:00 UT on February 15, approximately 2 hours before the X2.2-class hare was detected 
by HMI/FDO. It covers a 216 x 216 (Mm^) region, divided into a 600 x 600 grid. It is 


obtained using the very fast inversion of the Stokes vector (VFISV) algorithm (Borrero et ah 
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2011) based on the Milne—Eddington approximation. A minimum energy method (Metcalf 


1994 Leka et ah 2009) was used to resolve the 180° ambiguity in the azimuth angle of the 


magnetic held. In this study, the vector held is preprocessed in accordance with Wiegelmann 


et al. (2006). 


A numerical box with dimensions of 216 x 216 x 216 (Mm^) is given by 1 x 1 x 1 
as a non-dimensional value. 2x2 binning of the original data yields grid numbers which 
are assigned as 300 x 300 x 300. The detailed numerical scheme is available in Inoue et al. 


(2013 

) or 

Inoue et al. 

..Q 

r—1 

0 

CM 


2.2. NLFFF Extrapolation 


We hrst perform an NLFFF extrapolation to obtain the 3D structure of the magnetic 
held prior to the X2.2-class hare, and to understand its physical properties relating to sta¬ 
bility; we apply the MHD simulation as an initial condition. This is based on the MHD 


relaxation method developed by Inoue et al. (2014b) which shows the detailed algorithm. 


The initial condition is given by a potential held obtained from the Green function method 


(Sakurai 1982). At the boundaries, the three components of the magnetic held are hxed; 


in particular, the three observed components of the magnetic helds are set at the bottom 
boundary and their velocities are set as zero. The Neumann boundary condition is imposed 
on (j). A formulation of r] is given by 


Vnifff = 5?7o + Vi- 


\J X B\\v\ 


\B\ 


( 6 ) 


where 7^0=1.0 x 10“® and 771 = 1.0 x 10“^. The NLFFF is performed in the center area of 
AR11158 to avoid any ehects from the inconsistent force-free a, and outside it is hxed by 
the potential held (see Figure 1(a) in Inoue et aL||2014^ . This calculation corresponds to 
Run A in Table 1 as shown in Inoue et al. ( 2014a[ ). Figure [^a) shows a distribution of the 
magnetic twist and 3D held lines structure in the central area of AR11158. Magnetic twist 
is dehned by 

1 r J\ 


T 

-*■ r 


= ^ 

Air \B\ 


(7) 


where = J ■ B/\B\ and dl is a line element (Berger & Prior 2006). We found that the 


strong magnetic twist of the NLFFF is distributed in the range from a half-turn to one-turn 


twist. Inoue et al. (2014a) found that the NLFFF exists in a stable state not only against 


an ideal MHD instability, such as a kink or torus instability, but also against even a small 
perturbation. In other words, this NLFFF needs to be excited into a higher energy level to 
release the free magnetic energy required to produce the hare. 
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2.3. MHD Relaxation and Simulation 


We further perform the MHD relaxation using the NLFFF as an initial condition to 
excite it into higher energy level. We employ the zero-/? MHD equations which are the same 


as those used in the NLFFF calculation, except the velocity limiter is released (see Inoue et 


ah 2014b) and the resistivity formulation is replaced by the anomalous one as follows, 




j < jci 

J ^ jci 


( 8 ) 


where 772 = 5.0 x 10“^, and jc is the threshold current, set to 30 in this study. The boundary 
condition is same as a case of the NLFFF. Because the anomalous resistivity can enhance a 
reconnection in the strong current region between the twisted lines formed in the NLFFF, 
we can expect it to form stronger twisted lines. This calculation corresponds to Run C in 
Table 1 as shown in Inoue et ah ( 2014a| ). Figuresj^b) shows the distribution of the magnetic 
twist and 3D field lines structure after the MHD relaxation process, these are same format 
with Figure [^a). We found that the strongly twisted lines with more than a one-turn twist 
are formed after the MHD relaxation process, which are not shown in the original NLFFF, 
as shown in Figure [^a). The formation of these strongly twisted lines is clearly due to the 
anomalous resistivity, which is reminiscent of the tether-cutting reconnection. 

Finally, we perform the MHD simulation using the magnetic held shown in Figure [^b) 
obtained after the MHD relaxation process. At the boundaries, the tangential components 
of the magnetic held are released to interact with an induction equation consistently, while 
the normal component of the held remains constant. Note that in this case, the tangential 
components gradually go back toward the state of the potential held (or another less twisted 
state) because the NLFFF cannot maintain the force-free state completely. The velocity 
and 0 are treated in the same manner for all boundaries as were the NLFFF calculation 
and time-relaxed simulation. This calculation corresponds to Run D in Table 1 as shown in 
Inoue et ah (2014a). The magnetic held shown in Figure [^b) is already in a non-equilibrium 
state;therefore, it illustrates the dramatic dynamics shown in Figure [^c). An overview of 


these is was summarized in Inoue et al. (2014a). 


Results 


Inoue et al. (2014a) presented the overview of the dynamics and compared these with 


some observations. In this paper, we study the more detailed 3D dynamics, in particular, 
focusing on the formation of the large hux tube, he., the growing process of the CME, and 























other phenomena associated with the reconnection, eventnally nnderstanding a relationship 
flares and CMEs. 


3.1. Temporal Evolution of the Magnetic Twist 


In general, 3D dynamics obtained from the MHD simnlation show a complicated behav¬ 
ior, so it is not easy to extract essential components from them. However, the magnetic twist 
deflned in eqnation Q is one solid tool that we can nse to extract these components, e.g., a 
dynamics of flnx tnbe and the reconnection associated with the formation and dynamics of 
the flnx tube. Because the flux tube consists of a bundle of fleld lines with a strong magnetic 
twist, we can trace the temporal evolution of the flux tube by tracing that of the magnetic 
twist. Furthermore, since the magnetic twist is a magnetic helicity generated by a toroidal 
current in the flux tube (Berger & Prior 2006 Torok et ah 2010; Inoue et ah 2012 ), the 


magnetic twist would help us understand the reconnection dynamics in terms of conservation 
law of the magnetic helicity. 

We first show in Figure the temporal evolution of the magnetic twist during the flare 
obtained from the MHD simulation, mapped on the solar surface where T„ = 1.0 and T„=0.5 
are plotted. The locations surrounded by the contour of T„=1.0 and 0.5 correspond to where 
the base of the flux tube exists. We clearly see that the locations of the twisted flux tube 
are not flxed;rather, they move in a course. At positive polarity, the region occupied by the 
strong magnetic twist is moving southward. The region located at the negative polarity is 
shrinking in an early phase and then seems to maintain its size in a later phase. These are 
not simple pictures that we can image easily, particularly since they are occurring in 3D 
space. Furthermore, although the twist is widely distributed in the area close to the polarity 
inversion line (PIL) at initial time, it seems to be concentrated in the area away from the 
PIL at t=12.5. In other words, most of the strongly twisted lines disappear close to the PIL. 
The reason for this is that the post-flare loops are formed through reconnection during the 


flare, as shown in Figure 8(a) of Inoue et ah (2014a). 


3.2. Formation of the Large Flux Tube 

We see the temporal evolution of the magnetic fleld lines in Figure where we focus 
on the strongly twisted lines in red which have more than one-turn twist and those in blue 
which have less than a one-turn twist (since most parts are less than half-turn, hereafter we 
call them weakly twisted lines). These weakly twisted lines are traced from the location close 
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to where the strongly twisted lines exit. Although the strongly twisted lines with more than 
a one-turn twist launch away from the solar surface, they split into two parts in an early 
phase, f=2.5; one footpoint of some field lines is rooted in the negative polarity at southeast 
area which locates away from the central area. That is why we see that the strongly twisted 
region occupied on central negative polarity shrinks (see Figure]^. Interestingly, the weakly 
twisted lines enter the large flux tube, while the short twisted lines marked by the arrow in 
Figure]^ are left where the strongly twisted lines existed at f=0.0. The contour of the strong 
twists (T„=1.0) is plotted on the solar surface at f=0.0 and f=12.5. At f=0.0 the footpoint 
of the red held lines is rooted in the region surrounded by the contour of T„=1.0, while at 
f=12.5 this contour surrounds the footpoint of the blue held lines. 


This result indicates that the twisted lines that were initially weak convert into the large 
hux tube having strongly twisted lines with more than a one-turn twist. This process might 
be explained by the tether-cutting reconnection between the eruptive strongly twisted lines 
in red and ambient weakly twisted lines in blue, because the weakly twisted lines can receive 
magnetic helicity from the strongly twisted lines via reconnection processes. Therefore, the 
strongly twisted region on the central positive polarity moves the southern region occupied 
by the weakly twisted lines as seen in Figure 

Figure also shows the temporal evolution of the magnetic held lines from another 
angle, plotted in the same format as Figure except that the evolution of the strongly 
twisted lines is omitted. The yellow surface corresponds to the critical height of the torus 
instability, he., a decay index of n=1.5 where n is dehned by 


n{z) 


z d\B 


(9) 


where, the decay index is calculated based on the potential held. Because the decay index is 


derived from the external poloidal held sustaining the hoop force by the hux tube(Torok & 


Klieml 2007), it is difficult to separate the external held from the NLFFF, following Fan & 


Gibson (2007) and Aulanier et al. (2010). This result clearly shows that the weakly twisted 


lines reform into large and strongly twisted lines through reconnection before reaching the 
critical height of the torus instability. After this, the top of the held lines exceed the critical 
height and would further grow into the CME. 


We carry out a more detailed quantitative analysis in terms of the conservation of 
magnetic helicity. Figure [^a) displays the distribution of the magnetic twist at f=0.0 with 
the contours of T„=1.0 at f=15.0 and Bz. Although one footpoint of the large hux tube 
formed during the hare is observed within a region surrounded by the red line in Figure 
l^a), the twist value of the held lines at this location at f=0.0 is less than 0.375 turns, 
approximately. The twisted held lines at f=0.0 are plotted in Figure |^b), where the red 
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twisted lines have more than a one-turn twist while one of the blue twisted lines, ’WLl’, has 
less than 0.375 turn twist approximately and another blue twisted line, ’WL2’ has a twist 
close to zero. However, very few field lines in ’WL2’ have more than a half-turn twist, as 
shown in Figure |^c). If the tether-cutting reconnection takes place only between the field 
lines ’WLl’ and ’WL2’ , then the long twisted lines with more than a one-turn twist are 
never produced, since by the conservation law of the magnetic helicity. WLl and WL2 do 
not supply the twist enough to form such lines. This large flux tube is, therefore, formed 
through the reconnection between the weakly twisted lines ’WLl’, and ’WL2’, and strongly 
twisted lines which are a source supplying the magnetic helicity. 

The field lines at t=0.0 and t=15.0 traced from the same location on the positive polarity 
are plotted with distribution and contours of T„=0.5 and 1.0 at f=15.0 in Figures |^d) 
and (e), respectively. As we presented in Figureat f=15.0 the blue twisted lines become 
strongly twisted lines, even though they were weakly twisted at t=0.0. Eventually, the 
footpoint of the red held lines dominates at f=15.0 at a location marked by the dashed 
white circle, which is where the blue held lines resided at f=0.0. These results also support 
that magnetic reconnection takes place between the strongly twisted lines and weakly twisted 
lines, WLl and WL2;consequently, part of the red lines remain at the solar surface at f=15.0 
as sheared post-hare loops. Thus, some of the magnetic helicity accumulated into the strongly 
twisted lines at f=0.0 is transferred into the weakly twisted lines through the reconnection, 
resulting in the formation of the large strongly twisted lines. 


3.3. Reconnection Dynamics in the X2.2 Solar Flare 


We investigate the reconnection process in the early hare phase. Following [Toriumi 


et ah (2013) and Inoue et ah (2014a), we estimate the reconnected held lines by using a 
spatial variance of the held line connectivity by allowing them only by reconnection. The 
formulation is dehned by the following equation: 


d(XQ, tn) 1^1 ^n+l) ^l(^05^n)| (^n = 0.3), 


where Xi(xo, tn) is the location of one footpoint of each held line at time which is traced 
from another footpoint at Xq. Eventually, we calculate 


A(xo,t) 


6 (iUO) 


( 10 ) 


where the enhanced region in A{xo,t) indicates memories in which the reconnection took 
place dramatically in the twisted lines. Figure [^a) shows the distribution of A{xo,t = 2.0) 
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over the -B^(>0.01) distribution, which we previously discussed, and this shape is similar to 
that of an observed two-ribbon flare ( Inoue et al.||2014a[ ) Figure [^b) shows the A map with 
the twist contours T„=0.5 and 1.0 over the Figure [^a). From this hgure we found that the 
enhanced region ’T2’ is found just outside of the region containing the strongly twisted lines, 
while the region ’Tl’ overlaps with their insides. We can hnd the same prohles on the positive 
polarities. Figure |^c) shows the 3D magnetic held lines, where the red held lines indicate 
the strongly twisted lines with more than a one-turn twist and where the orange lines go 
through an inner part of the region surrounded by the current density | J|=5.0 contour. In 
an earlier hare phase the strongly twisted lines with more than a one-turn twist erupt away 
from the solar surface, beneath which the orange held lines are reconnected, forming the 
long twisted lines and post-hare loops on the A map. 


We investigate the relationship between the held lines structure and EUV images taken 
by SDO to further understand the reconnection dynamics. Figures [^a)-(c) represent the 
held lines structure at t=0, where each colored held line( red, orange, and blue) is traced 
from the regions, Rl, R2, and R3 surrounded by dashed circles in Figure [^a). The intensity 
contour taken from an AIA 171 A image is also plotted with these held line structures, where 
Figures [^a) and (b) show the much earlier hare phase and Figure [^d) exhibits it in middle 
hare phase. In Figures [^a) and (b), the most enhanced areas are located at El and E2 for 
each time where the small hux tubes are coexisting, e.g., we can see the red and orange held 
lines at El and three diherent held lines at E2. This means that the reconnection might 
take place among these small hux tubes. Figure [^d) shows the held lines structure t=2.0. 
The held lines traced from the regions Rl and R2 convert into post-hare loops, on the other 
hand, some of footpoint in blue traced from R3 are rooted into the hook region on another 
polarity in which the contours of the strongly twisted lines are also observed, as shown in 
Figure l^b). Another footpoint is rooted into the negative polarity located southeast away 
from the central area. This location is marked by a solid white circle at each time in Figures 
[^a)-(c), in which EUV images are enhanced during the hare. These results demonstrate 
that the location close to the PIL is crowded with some twisted lines, i.e., small hux tubes, 
before the hare, and then the complicated reconnection takes place among them during the 
hare. 


We also checked the A map and 3D magnetic held at t=4.0. Figure [^a) shows the 
A{xo,t = 4.0) with the twist contours mapped on the bottom surface in the same format 
as that in Figure [^b), and it well resembles the two-ribbon hare structure (see |Inoue et al.| 
2014a[ ). Figure l^b) exhibits the 3D structure of held lines in the same format as Figure [^c), 
except the distribution of the vertical velocity is plotted on the vertical cross section. Figure 
[^c) shows the enlarged view of the current layer and the distribution of vertical velocity 
close to the solar surface. We clearly see the upward and downward hows, which have an 
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origin at the inner side of the current layer, meaning that the reconnection is occurring here. 
The strongly twisted lines in red are a slightly different from that seen in Figure [^c). We 
can see that one footpoint is rooted into the negative polarity located southeast away from 
the central region. This result indicates that the strongly twisted lines also interact with 
the ambient twisted lines during their ascension, as suggested by Figure The orange held 
lines show similar behavior as that seen in Figure |^c). Through reconnection, the held 
lines convert into the long twisted held line, whose one footpoint is rooted into the negative 
polarity in the southeast, and the post-hare loops. Figure [^b) also clearly shows that the 
locations of both footpoints of the post-hare loops correspond to those in which the A maps 
exist. 


3.4. Comparisons with the Observations 
3 . 4 . 1 . Enhancement of the Tangential Field at Solar Surface 


In this section, we compare our simulation results with some observations to evaluate 
the validity of our reconnection dynamics. 


Wang et al. (2012) discovered an enhancement of the tangential magnetic helds during 


the hare, residing in a location close to the PIL on the solar surface. They, as well as Sun et 


al. (2012), suggested that this enhancement is strongly related to the reconnection. We check 


this with the temporal evolution of the tangential helds + B^) obtained from our 

simulation. Figure [^a) shows the B^ distribution in the central area in which the tangential 
helds are measured. Figures [^b) and (c) represent the distribution of the tangential magnetic 
helds at f=0.5, and f=15.0, respectively. We clearly see that the tangential helds are strongly 
enhanced at a location close to the PIL, whose value can reach more than 2000 G. These 
values and locations are therefore well consistent with those observed by Wang et al. (2012). 


The structure of held lines is depicted in an early hare phase above the enhancement 
of Bt at the photosphere, at f=0.0, 1.5, and 3.0 in Figures |^d), respectively. At f=0.0, just 
after the MHD relaxation process, the long sheared held lines are formed, and then they 
convert into the short loops at t=1.5. Since these short held lines pass through the strong 
current region, we can interpret that these correspond to post-hare loops formed through 
the reconnection. Furthermore, Bt is enhanced under the area in which the post-hare loops 
reside, even in the early phase at t=1.5. At t=3.0, we found that Bt is more enhanced 
and the taller post-hare loops are formed, above which the further new post-hare loops are 
being produced and piled up through sequential reconnection. This result implies that the 
tangential held of the newly formed post-hare loops sufficiently compresses the pre-existing 
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ones below, resulting in the enhancement of the tangential magnetic helds near the solar 
surface. 


In Figure [^e), we show a temporal evolution of an averaged shear angle between the 
tangential fields in the simulation and that of the potential field, which is measured on 
Bt > 0.4 at the photosphere, surrounded by white dashed square in Figure |^b). The shear 
angle is defined by V’ = cos~^{BfBpt)/\Bt\\Bpt\ where Bt = By), and Bpt = {Bp^, Bpy) 
(tangential components of the potential held). Following this result, the averaged shear angle 
< -0 > is decreasing, i.e., relaxing toward the potential held. These results indicate that the 
magnetic helds become overall less energetic as shown in Su et al. (2007) even though Bt is 
enhanced with time. 


3.4-2. Enhancement of the normal current density at the Solar Surface 


Janvier et al. (2014) reported an enhancement of the normal current density (J^) on 


the solar surface during the hare, particularly at the negative polarity in the central area 
(Region S- shown in Figure 3 in their paper). We also check the temporal evolution of the 
negative current density hux to compare with their observational result measured in the 
region S-. Figure [l0|(a) shows a distribution of at f=0.5 and orange line corresponds to a 
contour of Jz=-35 at which the Jz is strongly enhanced. The normal current hux is measured 
in an area surrounded by red dashed square in Figure [l0|(a) which almost touches to the 
enhanced region of Jz and close to the region S- set in Janvier et al. ( 2014[ ). A temporal 
evolution of the negative normal current hux J | — Jz\dS is plotted in Figure [To|(b). Our 
simulation results also shows its sudden enhancement, but its value decreases gradually after. 
This is diherent from the results found by [Janvier et ah' (2014) who observe a saturation 
of the normal current density. Figure [Io|(c) shows the temporal evolution of the contour 
for Jz = —35. This contour is moving away from the PIL as time progresses as shown in 


Aulanier et al. (2012), and the size becomes smaller as speculated in Figure dOla). 


We also checked the relationship between the enhanced current density and held line 
structure. The structure of held lines at t=3.5 passing through the inside of the current 
layer |J| = 10 is plotted in Figure [T^d). The one footpoint of these lines is rooted in the 
negative current region, where the current density is strongly enhanced. An insert exhibits 
the top view of the Jz map at t=3.5 where the arrow points out the most enhanced region 
of the negative current density. We further calculate the norm to elucidate the property of 
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these field lines. The norm is defined by following eqnation eqnation (Demonlin et al.|1996 ) 


N{x,y) = 


i 


E 

i=l,2 


dX, 

dx 


+ 


dX, 

dy 


( 11 ) 


where (Xi, X 2 ) is the relative distance corresponding to (x' — x", y' — y"). (x', y') and (x", y") 
are the positions of the end points of the field lines whose starting points are two adjacent 
grid points located at (xg, y'o) and (xg, j/g) on the photospheric snrface. This means that the 
locations of the end points of these field lines, which are traced from these starting points 
across a large N(x,y) valne, may differ greatly. Fignre [T^e) plots the field lines at t=3.5. 
These are the same field lines as in Figure [T^d), with a distribution of N(x,y) added. We 
can clearly see that the one footpoint of these field lines is rooting into the enhanced negative 
current region and is across the N(x,y) marked by Ql. This means that the footpoints of the 
field lines in a different topology defined by N (x,y) coexist in the enhanced current layer. We 
further plot the field lines at t=5.0 in Figure [^f), when all of field lines traced from the same 
location in Figure [T^e) become post-flare loops. From these results, the field liens plotted in 
Figure [lOl(d) or Figure [To|^e ) are in a sate just before and after the reconnection, and we found 
that the normal current density is strongly enhanced at the moment of reconnection, where 


the footpoint is rooted. These results are consistent with Janvier et ah (2013) demonstrating 


a relationship between the quasi separatrix layer (QSL) and current ribbons formed at the 
photosphere reproduced in their simulation. 

In our simulation, the normal current flux cannot maintain the peaking value in its 


temporal evolution. This is in contrast to the observation made by Janvier et ah (2014). 


One possible reason for this difference might be related to our boundary condition arising 
from the initial NLFFF. The temporal evolution of the tangential fields in our simulation 
relaxes the shear gradually through the photosphere such that it trends toward the potential 
field or another low energy state gradually so that cannot trace the evolution as seen in the 
observation. Also the gap between the observation and simulation shown in Figure [^e) might 
be caused by this problem. On the other hand, following Janvier et al.| ( |2013[ ), the current 
ribbons are linked with the current region in the solar corona at which the reconnection 
taking place, such that numerical diffusion might weaken this link. 
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4. Discussion 

4.1. Change for the Connectivity in the Strongly Twisted Lines 

We found that one fooptpoint of the new large flux tube that had the strongly twisted 
lines and T„ >1.0 and was formed through reconnection during the flare was rooted into the 
negative polarity located southeast away from the central area as seen in Figures [^d) or 
Bb). We further check it by comparing it with another AIA images. Figure [T^a) exhibits 
the AIA 171 A taken by the SDO at 01:48:25 UT on 2011 February 15, corresponding to the 
onset of the solar flare. The region surrounded by the black square corresponds to the area 
in which one footpoint of the flux tube was rooted during the flare, while another is rooted 
into the central area. The temporal evolutions of this area are shown in Figures [IT|b)-(e). 
As time progresses during the flare, we can see that the dimming region is growing. Figure 
[I^a) shows another AIA image, taken during the middle flare phase and superimposed with 
Bz contours. The black dashed circle corresponds to the region surrounded by the square 
marked in Figure [IT|a). The strongly twisted held lines at t=15.0 are plotted in Figure 

at t=15.0 in white on the positive polarity. We can see that another footpoint is rooted into 
the dimming region at the negative polarity, which is marked by the dashed circle in Figure 

Hla)- 

Figures [T^c) and (d) map the distribution of the open-closed held lines obtained from 
the NLFFF before and after the flare, respectively, with the held lines traced by the blue 
solid circle. At the region marked by the blue circle in both panels, the open magnetic held 
lines going through the side or top boundaries of this held of view (FOV) dominate before 
the flare, while the closed held lines with both footpoints are anchored within this FOV 
dominate after the flare. Since this process implies the reconnection between the open and 
closed held lines and produces the large flux tube as shown in Figure]^ this might support 
the dynamics obtained from our MHD simulation. 


12 b), where these held lines are traced from the inside of the contour for the twist T„=1.0 


4.2. A Relationship between 3D Dynamics and Enhancement of the 

Tangential Fields at the Solar Surface 

In section 3.4.1, we saw an enhancement of the tangential magnetic flelds at the solar 
surface that reside close to the PIT. In this section, we explore a relationship between this 
enhancement and the structure seen in the 3D volume. We first plot the temporal evolution 
of the maximum value R and summation (R (=WR)), for which the Bt >0.8 is 
counted, of the tangential magnetic flelds in Figure [I^ a). These values first increase during 
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the initial flare phase (from t=0 to t=5.5), where the increasing ratio is 38% measured in 
but these values saturate approximately in the late phase, as seen by IWang et ah 


(2012). This transition might be explained by different 3D dynamics in the early and late 
flare phases. We therefore check the 3D magnetic structure to explain this transition. 


Figure [I^b) shows the 3D view of the magnetic held at f=3.5, 5.0, and 6.5, respectively, 
with current density contours \J\ = 5.0 and isosurfaces corresponding to the critical height 
of the torus instability. At f=3.5, when and are increasing, the weakly twisted 

lines in blue which become the strongly twisted lines at f=6.5 are passing through the inside 
of the strong current region. This means that the large flux tube is forming in this phase 
through reconnection. At f=5.0 the large flux tube has been formed and its top is close to 
the location of the threshold of the torus instability. The profile of Bt^^^ or Bt^^^ are just 
about to begin their saturation phase. Eventually, the large flux tube exceeds the threshold 
of the torus instabilityflhen, Bt^^^ or Bt^^^ remain in the saturation phase. These results 
imply the profile of Bt shows two different picture in the formation of the large flux tube and 
its ascending phase over the threshold of the torus instability. The different 3D dynamics in 
the solar corona might be shown at the solar surface in the form of the enhancement of Bt. 


Essentially, this would be controlled by the reconnection in the solar corona. In the ear¬ 
lier flare phase, the reconnection takes place in the lower corona, meaning that the magnetic 
held is reconnected effectively and can enhance the tangential magnetic held at the photo¬ 
sphere because the strong magnetic flux density exists there. On the other hand, at a later 
phase, the magnetic reconnection point shifts toward the higher coronal region associated 
with the rise of the flux tube, so the reconnection is not working effectively relative to the 
one in the lower corona, which might put the brake on the enhancement of the tangential 
fields at the photosphere. These more detailed and quantitative analyses not presented here 
are left as future work. 


4.3. A Transition from Flare to CME 


Figure lA shows a summary of the dynamics for the X2.2-class solar flare in AR11158 
based on our simulation results. The onset of the solar flare starts in phase (i) where the 
strongly twisted lines are produced via tether-cutting reconnection of the twisted lines formed 
in the NLFFF. Consequently, the strongly twisted lines erupt away from the solar surface, 
and the sheared two-ribbon flares are appearing and new large flux tube is being formed in 
phase (ii). The CME onset starts in phase (iii), where the newly formed large flux tube is 
ascending upward over the threshold of the torus instability. 
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Inoue & Kusano (2006) suggested an existence of a critical height above which a flux 
tube cannot recover to the equilibrium state. They concluded that the CME formation 
due to the eruption in the flux tube is not enough to surpass the critical twist due to the 
kink instability;however, it is important to couple this with the loss of equilibrium such 
that the flux tube can grow over a certain height. Reaching the critical height of the torus 
instability would equate to the point where equilibrium is lost. In this study, the tether¬ 
cutting reconnection plays an important role as driver in the lower corona at the initial 
phase of flare-CME dynamics, rather than the kink instability. However, a scenario of the 
flare-CME relationship is basically the same as in our previous study Inoue & Kusano (2006). 


5. Summary 

We investigated the dynamics of the magnetic field during the X2.2-class solar flare 
taking place in AR11158, in particular focusing on the large flux tube formation and the 
related phenomena associated with the flares, eventually suggesting a relationship between 
flares and CMEs. 


1. We first reconstruct the NLFFF based on the photospheric magnetic field prior to the 
X2.2-class solar flare in AR11158 on 2011 February 15. This NLFFF never showed 
the dramatic dynamics as seen in observations because it is stable even against small 
perturbations (see Inoue et al.j2014a ), he., it is stable against the ideal MHD instability 
such as like kink and torus instabilities. The NLFFF, therefore, needs to be excited 
into higher energy levels to cause the huge flares. 


2. We further performed the MHD relaxation, he., the MHD simulation with the hxed 
boundary condition where the three components of the magnetic held are hxed at the 
bottom surface, using the NLFFF as an initial condition. An anomalous resistivity 
is applied in this calculation and as a result, the strongly twisted lines are produced 
through the reconnection between twisted lines taking place in the strong current 
region formed in the NLFFF. Because this magnetic conhguration cannot keep the 
equilibrium state, the strongly twisted lines successfully showed the eruption due to a 
loss of equilibrium. In fact, the same kinds of trigger processes are needed, rather than 
the anomalous resistivity, to produce the strongly twisted lines through the tether¬ 
cutting reconnection and bring them up to a critical height of the torus instability. 
We here point out some pioneering works about the triggering process of the solar 


flare: Flux cancellation on the photosphere ( 

van Ballegooijen & Martens 

1989 

Green 

et al.|2011 

Amari et al.pOlO; 

Amari et al.||2011 

), emerging magnetic flux through the 

photosphere ( 

Heyvaerts et al. 

1977 

Feynman & Martin 

1995 

Chen & Shibata 

2000 
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Roussev et al. 

2012 

Kusano et al.||2012 

Toriumi et al. 

2013 

Bamba et al.|2013 


processes should be taken in account in our simulation as important future work so 
that our model does not depend on the anomalous resistivity model. 

3. Although the strongly twisted lines with more than a one-turn twist launch away 
from the solar surface, they further reconnected with the ambient weakly twisted lines. 
Consequently, the new large flux tube is formed. This result shows that the large flux 
tube growing into the CME is formed through the tether-cutting reconnection between 
the strongly and weakly twisted lines. The tether-cutting reconnection, therefore, plays 
an important role in the triggering process of the CME, he., the formation process of 
the large flux tube and overcoming the threshold of the torus instability, as well as the 
triggering process of the flare, he., the formation of the strongly twisted lines in the 
NLFFF driving the initial eruption. Our results suggested that the flare and CME are 
connected via these two scenarios. 


4. We investigated a relationship between the distribution of A dehned in equation (10), 
corresponding to the connectivity change of a twisted held line, and the distribution of 
the twisted lines. This A distribution well resembles the distribution of the observed 
two-ribbon hare in the earlier hare phase. Most of the strong A region is found just 
outside of the region of the strongly twisted magnetic held. On the other hand, we 
found that the footpoint of the long twisted lines is anchored into the hook area formed 
in the A distribution surrounding the footpoint of the eruptive strongly twisted lines. 
This scenario would be similar to the 3D hare model in Aulanier et ah (2012) and 


Janvier et al. (2013). However, following this point, the eruptive strongly twisted lines 


showed complicated behaviors, in particular, the reconnection between the strongly 
and weakly twisted lines resulting from the complex magnetic topology of AR 11158, 
through which the large hux tube is formed as mentioned above. 


5. In our simulation, we observed a portion of the tangential magnetic held, Bt at the 
solar surface residing close to the PIL in the earlier hare phase, whose maximum value 
was recorded to be more than 2000 G. During this phase, the dynamics of the twisted 
lines are observed within a region less than the critical height of the torus instability in 
which the large hux tube formation is beginning. We further found that the post-hare 
loops already exist and cover the region on the photosphere where Bt is enhanced, even 
though the shear angle between the tangential components in the simulation and that 
of the potential helds decreases there gradually with time. Therefore, this enhancement 
is inferred due to be piled up by the new post-hare loops produced by the sequential 
reconnection, he., Bt on the photosphere is compressed by them and then enhanced. 
After the large hux tube forms and exceeds a critical height of the torus instability. 
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the value of approximately reaches saturation. The maximum value increases 


during the initial flare phase while Wang et ah (2012) reports the average value at the 
selected region increases by ~ 30%. These results indicate the different prohle of Bt on 
the photosphere, such that increasing or saturating behaviors might be projected by 
the different 3D dynamics observed in the solar corona. On the other hand, although 


Wang et ah (2012) and Sun et ah (2012) pointed out the tangential helds become better 


aligned with the PIT during the flare, these ’more parallel held lines’ appearing in the 
lowest area were not reproduced in our simulation. 

6 . We also found the enhancement of on the solar surface in our MHD simulation. 


which is consistent with the analysis of the observed result by Janvier et ah (2014). 


This enhancement is observed in a region where a footpoint of the held lines in which 
reconnection is taking place in the strong current region is rooted, as pointed out by 


Janvier et ah (2013). On the other hand, our MHD simulation showed that the value of 


Jz is decreasing gradually as time passes after its enhancement, which is diherent from 
the observed picture. There are still some gaps which remain between the observations 
and simulation. 


In this paper we investigate the dynamics of magnetic helds during the X2.2-class hare 
taking place in AR11158, focusing in particular on the transition dynamics from the solar 
hare to the CME. Unfortunately, we cannot trace the evolution of the hux tube for a long 
enough time within the limited volume of the simulation, he., we cannot conclude whether 


or not the large hux tube shown in our MHD simulation is in its hnal form. Hu et ah (2014) 


reported a very highly twist hux tube traveling an interplanetary, up to 5 turns per AU, 
which decreases toward the edge. At least, our NLFFF or hux tube just launching from 
the solar surface whose twist values are concentrated mostly in a range from one turn to 


1.5 turns cannot explain this highly twist. Yamamoto et ah (2010) also made an ehort to 
clear the gap of twist between measured at the solar surface and interplanetary. But the 
problem has been left yet. A large-scale simulation covering the interplanetary space would 
be important to hll this gap. 

On the other hand, the twist accumulated in the large hux tube looks weak compared 
to the initial twisted lines in red in Figure |^a) even though the hux tube exceeds a one- 
turn twist. These would be caused by two problems: one is due to a numerical dihusion 


(Su et ah 2011); another is due to the initial NLFFF. In order to settle these problems, a 


high-resolution simulation is important to suppress the numerical dihusion, he., suppress the 
dihusion of the twist accumulated in the hux tube. Furthermore, a higher accurate NLFFF, 
achieving more close to an equilibrium state, would be also required to avoid release of a shear 
through the boundary. These problems must be settled to complete the data-constrained 
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or -driven simulation for solar flares and CMEs. On the other hand, future observations 
plan to observe a chromospheric magnetic held {e.g., SOLAR-C mission Watanabe 2014). 
Because a chromospheric magnetic held is in a lower f3 state than that of the photospheric 
magnetic held, it would help us dehnitely construct higher accurate NLFFFs and perform 
data-constrained or -driven simulations. 
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Fig. 1. — (a) The left panel is a distribution of the magnetic twist of the NLFFF 2 hours 
approximately before the X2.2-class flare, with contours of where black and white contours 
correspond to 5^=0.25(625 G), and -0.25, respectively. The held lines begin in the region 
(0.35Lo <x< 0.65Lo and 0.4Lo <y< O.TLq, where Lq is 216 Mm ). Then the magnetic twist 
is mapped, where the value of 15^1 is more than 0.01 (=25 G). The middle panel shows the 
held line structure. The orange curves show twisted lines with more than a half-turn twist 
while the blue lines have less than a half-turn twist. The Right panel also represents the 3D 
structure of the held lines from another angle, (b) The twist distribution and 3D held line 
structure after the MHD relaxation process. The hgure format is the same as (a), except the 
red lines show the strongly twisted lines with more than a one-turn twist, (c) The temporal 
evolution of the magnetic held obtained from the MHD simulation. The orange and blue 
lines have twisted lines more and less than a half-turn twist, respectively, at an initial time 
f=0.0. 
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- 0.25 „ 0.25 


Fig. 2.— Temporal evolution of the magnetic twist mapped on the solar surface, focusing 
on the central region of AR11158. The distribution is represented in gray scale. The 
red and orange lines are contours of the magnetic twist = 1.0 and = 0.5, respectively. 
Note that the inside regions surrounded by the red and orange lines are occupied by the 
strongly twisted lines. 
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Fig. 3.— Temporal evolution of the dynamics of the selected magnetic held lines with 
distribution. The red and blue lines represent the held lines more and less than a one-turn 
twist, respectively, at an initial time t=0.0. Note that most of the blue lines have weak 
twist less than 0.35 turns. The red arrow indicates sheared-post hare loops. The contour of 
r„=1.0 is plotted at surface in red at f=0.0 and 12.5. 
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t=0.0 (b) 


t=2.5 
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t=5.0 (d) 
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Bz 

Fig. 4.— Temporal evolution of the dynamics of the magnetic held lines from side view 
of Figure The held lines format is same in Figure except the temporal evolution of 
the initial strongly twisted lines with more than a one-turn twist is not plotted. The yellow 
surface corresponds to a critical height of the torus instability, i. e, the isosurface with decay 
index ?7,=1.5 where n is dehned by equation (§. Arrow points towards the intersection of 
the magnetic held lines and the threshold of the instability. 
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(d) (e) 
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Fig. 5.— (a) The distribution of the twist at t=0.0 with contours of T„=1.0 at t=15.0 in 
red and \Bz\=Q.2b in white, (b) The held lines at t=0 are plotted over (a). The format of 
these held lines are same as Figure]^ (c) The same held lines at t=0.0 are plotted over (a) 
except the strongly twisted lines in red are not plotted, (d) The twisted lines at t=0.0 are 
plotted in red and blue, traced from the positive polarity. The contours of twist T„=0.5, and 
1.0 at t=15.0 are plotted in red and orange, respectively, over the distribution where the 
twisted lines in blue are started from the inside region surrounded by the contour T„=1.0 at 
t=15.0. (e) The same format with (d) except the held lines are plotted at t=15.0. 
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(C) 



Fig. 6.— (a) A(a;o 5 2.0) is mapped on the distribution where A(a;o,2.0) is measured at 
\Bz\ > 0.01 and AO.3. (b) The contours of the twist T„ = 0.5 and 1.0 are plotted in the 
orange and red lines with A map over the B^ distribution, (c) The 3D magnetic held line 
structure is plotted over (b). The red lines represent the strongly twisted lines with more 
than a one-turn twist, while the orange held lines are passing through the region surrounded 
by the blue line, which is a contour of the current density (| J|=5.0). 
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(c) 


(d) 




-0.25 02 0-25 

Fig. 7.— (a)-(c) The field lines at t=0 are plotted with the contour of AIA intensity and 
A(a;o, 2.0) over the distribution. Each colored field line is traced from each colored region, 
Rl, R2, and R3 marked in Figure |^a). The purple lines correspond to the AIA intensity 
/o5f(/)=9.0 at 01:50:25 UT, 01:52:49 UT and 02:04:49 UT from (a) to (c) where the first two 
correspond to the much earlier flare phase and the latter is in the middle flare phase, (d) 
Each colored field line at t=2.0 is plotted over (a), but AIA intensity is not plotted. 
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Fig. 8.— (a) A(a;o 5 4.0) is mapped on the distribution at the solar surface with the twist 
contours, Tn =0.5 and T„=1.0 in orange and red, respectively, (b) The held lines structure 
at t =4.0. The hgure format is the same as that in Figure |^c), except the contour of the 
current density | J|=5.0 is plotted in white, and the vertical velocity distribution is plotted 
on the vertical cross section, (c) The vertical velocity close to the strong current region is 
plotted with the current density contour | J|=5.0. 
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Fig. 9.— (a) distribution with a polarity inversion line in white is plotted in the range 
of 0.33<x<0.63 and 0.37<y<0.67. (b)-(c) The distribution of the tangential components of 
the magnetic held Bt=-sjB1 + at f=0.5 and f=15.0. (d) The structure of held lines in the 
early hare phase is plotted in gray. The held lines shown in (e) and (f) are plotted passing 
through the inside of the current contour | J|=10. The held lines in (d) are traced from same 
points in (e). The bottom surface represents the distribution of which is in the same 
format as (b) and (c). (e) A temporal evolution of an averaged shear angle between the 
tangential helds in the simulation and that of the potential helds, measured on Bt > 0.4 at 
the solar surface, surrounded by white dashed square(0.433<x<0.533 and 0.50<y<0.567). 
depicted in (a). The unit is degree. 
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Fig. 10.— (a) A distribution of with contour of |i?2|=0.25 in white and Jz=—35 in 
orange at f=0.5. (b) Temporal evolution of the negative normal current flux, dehned by 
/I — Jz\dS, calculated in a region surrounded by red dashed square (0.433<x<0.49 and 
0.533<y<0.567) marked in (a), (c) Temporal evolution of the contour for the 7^=—35, at 
t=2.5(in red), 5.0(in green), and 7.5(in aqua), respectively, on distribution, (d) The held 
lines are plotted in red, passing through the inside region surrounded by the contour of | J| = 
10 over Jz distributions. Insert is the top view of the map where the arrow points out the 
most enhanced region, (e)-(f) The held lines are plotted in the same format as (c) and are 
plotted with a distribution A^(x, y) dehned in equation ( [II| ) at f=3.5 and t=5.0, respectively. 
The white contour represents the |i?2|=0.25. 
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Fig. 11.— (a) EUV image taken from AIA 171 A at 01:48:25 UT in the range of 216 x 
216 Mm^, corresponding to an area in which the NLFFF or MHD simulation is performed. 
(b)-(e) EUV images taken from AIA 171 A from t= 01:48:25 to t=02:ll:13 during the flare, 
in the range of a region surrounded by the black square in (a). 
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Fig. 12.— (a) EUV image taken by AIA 171 A in the range of 0.25<x<0.78, 0.33<y<0.67 
at 02:04:49UT. Red and blue lines correspond to the contours of R^=0.25, and -0.25, respec¬ 
tively. (b) Field lines with more than a one-turn twist at t=15.0, he., flux tube, are plotted 
over (a). White line corresponds to the contour of T„=1.0 at t=15.0 where the held lines 
are started from the inside of this region, (c) Map of open-closed held lines for the NLFFF 
on February 15 00:00UT corresponding to the 2 hours before the hare. White regions are 
dominated by the closed held lines where both footpoints are anchored in this held of view 
(FOV) (0.38<x<0.62, 0.4<y<0.63), while the black regions are dominated by the open mag¬ 
netic held whose one footpoint goes through the boundary of this FOV. The blue lines are 
held lines traced from the solid circle, (d) Map of open-closed held lines for the NLFFF on 
February 15 03:00UT, corresponding to an hour after the hare. Figure format is same as in 
(c). 
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Fig. 13.— (a) A temporal evolution of the maximum value of the tangential components 
of the magnetic held, and its summation of that, (=X] ^t) where Bt is counted 

over 0.8(=2000G), measured in the range of (0.25Lo <x< 0.75Lo, 0.25Lo <y< 0.75Lo) are 
plotted in red and blue, respectively. The vertical solid line indicates a location where the 
prohle of the Bf^^^ suddenly changes from an increasing phase to the saturating phase, (b) 
Temporal evolution of the magnetic held lines with the contour of \ J\ = 5.0 in white and the 
isosurface of the threshold of the torus instability. Bottom surface represents the distribution 
of Bt- The blue held lines are drawn in the same format as in Figure]^ and orange held lines 
pass through the region surrounded by the contour of the current density in white. Arrow 
pointed out the intersection between the magnetic held lines and isosurface. 
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Fig. 14.— Summary of the dynamics in AR11158, based on our NLFFF extrapolations and 
MHD simulations. The upper two insets show the 3D structure of the held lines in the same 
format as in Figure [T^b), except that the vertical velocity is plotted on the vertical cross 
section, where its top corresponds to the approximate critical height of the torus instability. 


















































